Selected topics in e+e - collisions 



Victor L. Chernyak, 

Budker Institute of Nuclear Physics, 630090 Novosibirsk, Russia, 
O ■ E-mail: v.l.chernyak@inp.nsk.su 



Content 



1) The leading twist pion wave function 
lO \ 2) "Improved" QCD sum rules with non-local condensates 

3) Pion and kaon form factors and charmonium decays: 
fSj ■ theory vs experiment 

^ . 4) 7*77r° - form factor 

! 5) The new non-local axial anomaly 

6) Cross sections 77 — > 7r + 7r~, K + K~ , K$Ks 

\0 • 1. The leading twist pion wave function 

O 

a: 



The general formula for the leading term of the hadron form factor in 
QCD was first obtained in pQ and has the form : 



(P2, s 2 , A 2 | J\\pi, si, Ai) = Ci 2 (-=)' 1 2 ' + 3 , (1) 

where g = Q»2 — Pi), w m i n is the minimal number of elementary constituents in 
a given hadron, n min = 2 for mesons and n min = 3 for baryons, Si j2 and Ai ;2 
are hadron spins and helicities, the current helicity A = (A1+A2) = 0, ±1, and 
the coefficient C\2 is expressed through the integral over the wave functions 
of both hadrons. It is seen that the behavior is independent of hadron spins, 
but depends essentially on their helicities. 

The largest form factors occur only for Ai = A 2 = mesons and Ai = 
A 2 = ±1/2 baryons. For two mesons, for instance: 



FiM) = F[f d \Q 2 ) (l + 0{a s ) + 0(k 2 QCD /Q 2 )^ 



1 Talk given at the International Workshop "e + e collisions from (f> to J/ip" 
March 1, 2006, Novosibirsk, Russia. 
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where /i j2 are meson coupling constants, for instance, f n ~ 130 MeV, fx = 
160 MeV etc., and 0i i2 (x, /l) are their leading twist wave functions which 
determine distribution of two meson quarks in fractions x and 1 — x of the 
meson momentum, they are normalized as: J dx<pi(x,~p) = 1. 

The properties of the pion wave function 07-(x, //) were investigated in 
[2] using QCD sum rules |3]. It was obtained that a few lowest moments 
of 4> n (x,fi ~ 1 GeV) are significantly larger than for the reference ("asymp- 
totic") wave function <p^ asy \x) = 6x(l — x), and based on these results the 
model form was proposed: 4>^ z \x^ = 30x(l — x)(2x — l) 2 , at the normal- 
ization scale /i = 0.5 GeV. This wave function is much wider than <p( asy \x), 
and this increases greatly calculated values of various amplitudes with pions. 

Recently there appeared first reliable lattice data on the value of a^/i = 
1 GeV), which is the second Gegenbauer moment of (f> n (x, /i = 1 GeV) : 2 : 



4(/x = IGeV) 



0. for 



x 



0(asy) ( 



0.50 for 



x 



^ z (x,// = lGeV) 



(0.19 ± 0.05) from "improved" QCD sum rules with 
non — local condensates, A. Bakulevet. al. jl] 



0.38 ±0.23 



~0.11 

-0.06 



lattice, L. Del Debbio et. al. [3] 



0.364 ±0.126) lattice, M. Gockeler et. al. ® 



It is seen that the lattice data disfavor <h(x,n = 1 GeV) ~ 4><f sy \x) and 
clearly prefer wider wave function, but it is highly desirable to increase the 



2 The original lattice results for a-i (/i) are obtained for /i 2 = 7.1 GeV 2 in [H] and 
fi 2 = 5 GeV 2 in 6 . They are evolved to fi 2 = 1 GeV 2 at NLO (see e.g. 0). 

Besides, the results for (£ 2 )t- in [H] are really obtained for the "pions" with the masses 
550 MeV < fin < 1100 MeV, and extrapolated then to the chiral limit ^ — > 0. It looks 
strange at the first sight that the values of (£ 2 ),- stay in nearly intact in the whole 
interval < fx v < 1.1 GeV, while one can expect that the wave function of the heavy 
"pion" with fj,^ = 1.1 GeV will be noticeably narrower (i.e. has smaller (£ 2 )) than those 
of the real pion. 



2 



accuracy of these lattice calculation (and to decrease //„.). 

2. "Improved" QCD sum rules with non-local condensates. 

The original (and standard) approach [3] for obtaining QCD sum rules 
calculates the correlator of two local currents at small distances expanding it 
into a power series of local vacuum condensates of increasing dimension: 3 

(01^(^(0)10) = £(z 2 rC< n (0|O n (0)|0>. (3) 

n 

In practical applications this series of power corrections is terminated after 
first several terms, so that only a small number of phenomenological param- 
eters (0|O n |0), n < n Q determine the behavior of many different correlators. 
This standard approach was used in |2j to calculate a few lowest moments 
of the pion wave function: (i^ 2 ™)^ = Jq dx (j) 7T (x)(2x — l) 2n . It appeared that 
the values of these moments are larger significantly than those for (f)( asy \x). 
The most important power corrections in these sum rules originated from the 
quark condensate ~ (0|g(0)g(0)|0) 2 . 

The " improved" approach [0] [TU] jl] proposed not to expand a few lowest 
dimension non-local condensates, for instance (the gauge links are implied) 
$(^ 2 ) = (0|g(0)rg(*)|0) and (0|G^(0)G M1/ (z)|0), into a power series in z 2 , but 
to keep them as a whole non-local objects, while neglecting contributions of 
all other higher dimension non-local condensates. This is equivalent to keep- 
ing in QCD sum rules a definite subset of higher order power corrections while 
neglecting at the same time all other power corrections which are supposed 
to be small. This is the basic assumption underlying this "improvement". In 
other words, it was supposed that the numerically largest contributions to 
the coefficients C n in eq.(3) originate from expansion of a few lowest dimen- 
sion non-local condensates, while contributions to C n from higher dimension 
non-local condensates are small and can be neglected. Clearly, without this 
basic assumption the "improvement" has no much meaning as it is impossible 
to account for all multi- local condensates. But really, no one justifications of 
this basic assumption has been presented in |§] [TU] [ij. 

3 In practice, the Fourier transform of eq.(3) is usually calculated, supplied in addition 
by a special " Borelization" procedure which suppresses contributions of poorly known 
higher dimension terms. This is not of principle importance, but is a matter of technical 
convenience and improving the expected accuracy. On account of loop corrections the 
coefficients C n depend logarithmically on the scale z 2 (or q 2 ). 
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Moreover, within this approach one has to specify beforehand not a few 
numbers like (0|G 2 (0)|0), (0|g(0)g(0)|0), but anumber of functions describing 
those non-local condensates which are kept unexpanded. Really, nothing def- 
inite is known about these functions, except (at best) their values and some 
their first derivatives at the origin. So, in [H] JU|[I] definite model forms of 
these functions were used, which are arbitrary to a large extent. The uncer- 
tainties introduced to the answer by chosen models are poorly controlable. In 
principle, with such kind of " improvements" the whole approach nearly loses 
its meaning, because to find a few pure numbers (£, 2n )n one has to specify 
beforehand a number of poorly known functions. 

As for (£, 2n ) n , the main "improvement" of the standard sum rules was 
a replacement (0|g(0)g(0)|0) 2 — > (0\q(0)q(x)q(y)q(z)\0) , "factorized via the 
vacuum dominance hypothesis to the product of two simplest (0|g(0)g(z)|0) 
condensates" 0^2] 111- Clearly, such a "functional factorization" looks very 
doubtful in comparison with the standard "one number factorization" of 
(0|g(0)g(0)g(0)g(0)|0). Using this approach, it was obtained in the latest 
paper |j: a 2 (/i = 1 GeV) ~ 0.19, a 4 (/i = 1 GeV) ~ -0.13. 4 This cor- 
responds to the effectively narrow pion wave function, with Jq 1 dx<f) 7T (x, /i = 
1 GeV)/x = 3(1 + a 2 + a^) ~ 3(1.06), i.e. nearly the same value as for the 
asymptotic wave function with a 2 = = 0. 

As for the above described basic assumption of "the improved approach" 
to the QCD sum rules, we would like to point out that it can be checked 
explicitly, using such a correlator for which the answer is known. 

Let us consider the correlator of the axial A u (0) = ^(0)7^75^(0) and 
pseudoscalar P(z) = u(z)i^z,d(z) currents: 

!„(*) = (0\TA,(0)P(z)\0) = ^I(z 2 ), (4) 



00 

i(z 2 ) = j:(z 2 rcnZ n (o\o n (o)\o), o , c n = £ c;=E^ Born,i) ( 1+ ^) ; 

n=0 i=a,b,c... i 

where Z n = Z n (u 2 ~ l/z 2 ,fi 2 ) are the renormalization factors of the opera- 
tors O n , while f n = 0(a> s ) is due to loop corrections to the hard kernels. 

4 For a not very clear reason this differs significantly from the previous results obtained 
within the same approach in the second paper in [Hj : 02 (m = 1 GeV) ~ 0.35, 0,4, {ji = 
1 GeV) ~ 0.23. 
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Figure 1: 



The exact answer for this correlator is well known in the chiral limit 
mu,d = 0, so that calculating it in various approximations one can compare 
which one is really better, and this will be a clear check. So, let us forget for 
a time that we know the exact answer, and let us calculate this correlator 
using the "improved" and standard approaches. 

The exact analog of the above described basic assumption predicts here 
that, at each given n, the largest coefficients in eq.(4) originate from the ex- 
pansion of lowest dimension non-local quark condensate Q(z 2 ) = (0\q(0)q(z)\0) 
shown in fig. la. Decomposing it in powers of z 2 , this results in a tower of 
power corrections with "the largest coefficients" C^ Born ' a ) : 

oo 

= J2(z 2 rcP°^(l + f:)Z n (0\O n (0)\0)^ o1 (5) 

n=0 

where O n are the corresponding local operators, O a ~ qq, 0\ ~ quGq, etc. 

The contribution of the fig. lb to I{z 2 ) in eq.(4) is originally described 
by the three-local higher dimension condensate (0\q(0)G(x)q(z)\0) . Its ex- 
pansion produces finally a similar series in powers of z 2 which starts from 
(Ox) and have coefficients Cp orn ' b )(l + /*). Besides, the dia grams fig.lc,... 
(not shown explicitly in fig.l) with additional gluons emitted from the hard 
quark propagator in fig.l produce similar series, starting from higher dimen- 
sion condensates (O n ),n > 2, and with the coefficients C^ Born ' c )(l + etc. 
In the framework of the above basic assumption, there should be a clear 
numerical hierarchy: Cp° m > a ) > C£ Bora > b ) > C^ OTn ^ > C(Born, i+i) ( all 
^Bom.i) are p aram etrically 0(1)), so that one can retain only the largest 
terms cK Born > a ) an d safely neglect all others. 

Let us recall now that the exact answer for this correlator is very simple 
(the spectral density is saturated by the one pion contribution only), and is 
exhausted by the first term with n = in eq.(4). 

In other words, there are no corrections in powers of (z 2 ) n>0 in this corre- 
lator at all. The reason is of course that other contributions with coefficients 
Q(y£a) j n eq.(4) neglected in the "improved" approach, cancel exactly all (ex- 
cept for the first one) "the most important" coefficients C^f 1 in eq.(5). For 
instance, the first power correction ~ (Oi) from the fig. lb diagram cancels 
the second term n = 1 in eq.(5). And power corrections from other higher 
dimension multi-local condensates from next diagrams not shown in fig.l, to- 
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gether with next corrections from the fig. lb diagram, cancel exactly all next 
n > 2 "the most important" terms from the fig. la diagram in eq.(5). 

So, the basic assumption of the "improved" approach clearly fails: those 
power corrections which are claimed to be "less important" in comparison 
with "the most important corrections", appeared to be not small and even 
cancel completely here all "largest corrections". 5 

5 When taking the Fourier transform (FT) of eq.(4) the relative values of terms with 
different n change, as only singular in z 2 terms contribute. So, in this simple correlator, 
only two first terms (O ) and (Oi) will contribute to FT of eq.(4) in the lowest order Born 
approximation: Z n — 1, C l n — C^° rn ' ^ = const„. All terms with n > 1 contribute to FT 
only on account of loop corrections. So, the contributions of n > 2 terms into FT of eq.(4) 
will have smaller coefficients ~ a s C^ mn \ 

In essence, all these technical complications are the extraneous issues for a check of the 
basic assumption of " the improved approach" , which is essentially the assumed hierarchy 
1°$ — (C r \ +1 /C7^) < 1 of coefficients C l n at each given n (and not the relative values of 
C n at different n) in the expansion of any correlator (0\TJi(z) J2(0)|0). And this basic 
assumption can be most easily and clearly checked just in the co-ordinate representation 
where the whole tower of terms survives in the correlator in eq.(4) even in the Born 
approximation. In any case, r4 — ci Born ' l+1 ^/ci Born ' , independently of whether it is 
calculated in the space-time or momentum representations, and so whether r n is small or 
not is independent of the representation and is a real check of the above basic assumption. 
Therefore, the statements like: "all terms with n > 2 enter the FT of eq.(4) with smaller 
(~ a s ) coefficients and are invisible in the Born approximation, so that a check of whether 
r n>2 i s smai l or n ot can't be performed" , look as an attempt to avoid a real check. 

In general, what number of terms in the series in powers of z 2 will survive after FT and in 
the Born approximation, depends on the correlator considered. In other correlators more 
condensates (O n ) will contribute to FT, even in the Born approximation. For instance, in 
the sum rules for (£ } w in :2J all (O n ) with dimO„ < (6 + 4fc) will contribute. So, even 
for (£ 2 )tt there will be a sufficiently large number of various condensates, even in the Born 
approximation. 

In m the appropriate FT of the correlator {Q\T d{y)^f ^u{— y) P(z)\0) was consid- 
ered within the "improved" approach, which resulted in the sum rules for the pion wave 
function ^> ff (x) itself (i.e. for all moments (£, 2k )-7r)- Accounted were only contributions 
from analogs of the fig. la and fig. lb diagrams, while all other contributions were finally 
neglected, according to the basic assumption. (Let us note that, even in the Born approx- 
imation, all condensates (O n ) with dimO„ < (5 + 4fc) contribute to FT in these sum rules 
for the moment (S,^)^)- Some model form was used for the quark bi-local condensate 
Q(v 2 ) — (0|g(0)q(w)|0) to describe the analog of the fig. la diagram, while the three-local 
condensate from the analog of the fig. lb diagram was (in essence, arbitrarily) "simplified" 
and expressed through the same &(v 2 ). The wave function (f> v (x) obtained in this way 
appeared to be significantly narrower that 0( as y)(x) (see fig.5 in [IT]), with ~ -0.26 
and (£ 2 )tt — 0.11 (somewhere at the low normalization point ^ ~ 1 GeV). These numbers 
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On the other hand, calculating the correlator in eq.(4) in the standard 
approach (which is a direct QCD calculation accounting for all terms of given 
dimension) , one finds that the sum of corrections of given dimension is zero, 
as it should be. 

The conclusion is that the above described "improved approach" to QCD 
sum rules can easily give, in general, the misleading results. 

3. Pion and kaon form factors and charmonium decays: theory 
vs experiment 

The calculated values of Br(%j — > 7r + 7r~) show high sensitivity to the 
precise form of <p n (x, fi ~ 1 GeV) |2] (see [12] for a review) : 



Br(x 



7T 7T 



3-10 



— i 



for 



x 



0(asy)( 



X 



1 • 10 



-2 



for 



x 



tCZ, 



X, Ho) 



k (0.50 ± 0.06) • experiment US 



(6) 



Br( 



f ~ 1 • 10" 4 for 



X 2 -> 7T 7T 



X 



t asymp , 



X 



0.24 • 10~ 2 for 



(X) = CZ (x,/i o ) 



(0.18 ± 0.03) • 10- 2 experiment [H] 



(7) 



It is seen that 4>)^ z '(x, /x G ) predicts branchings in a reasonable agreement 
with the data, while these numbers for (j) n (x) ~ (p^ asy \x) are not ~ 20%, but 
~ 20 times smaller. 

The pion form factor, see eq.(2), also shows high sensitivity to the form 
of (f) n (x,fj,). Below are given some typical numbers for F^ ead \\Q 2 \ ~ Q 2 Q = 

disagree both with the recent lattice results 0: afO" = 1 GeV) = (0.364 ± 0.126) and 
(£ 2 )^=iGcV = (0-325 ± 0.043), and with previous results of the same author [Hj. From 
our point of view, all this only illustrates that playing with various "improvements" of 
sum rules and/or with some model forms for non-local condensates, one can obtain very 
different results for (£, 2k )-n- 
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10 GeV 2 ) 6 , together with the old and recent data. 

' (0.69 ±0.19), Q 2 = 9.8GeV 2 , C.J.Bebek et. al. gH 

~ 0.13 GeV 2 for </y(x) = 0( as y)(x), a s ~ 0.3 

' ~ 0.5 GeV 2 for ^(x) = cz (x, // ), a s ~ 0.4 [2 

(1.01 ± 0.11 ± 0.07) GeV 2 , s = -Q 2 = 13.48 GeV 2 , 
CLEO - 2005 [T5] 

Estimates show that higher order corrections to the leading term contribu- 
tion to F n (Q 2 ) can constitute up to tens per sent at Q 2 ~ 10 GeV 2 . Besides, 
one naturally expects that the asymptotic behavior is delayed in the time-like 
region s = —Q 2 > 0, in comparison with the space-like region s = —Q 2 < 0, 
so that, for instance, \F n (s = 10 GeV 2 )/F 7r (s = — 10GeV 2 )| > 1, similarly to 
the nucleon form factor, see fig. 2 . In any case, this problem with the large 
value of \sF w (s = 13.5 GeV 2 ) | ~ 1 GeV 2 measured by CLEO is much 
more severe for (f) n (x) ~ 0( asj/ )(x), in which case the leading term contribu- 
tion has to be increased ~ 7 — 8 times to compare with the data, than for 
(f) n (x) ~ (j)^ z '(x), where the difference is about a factor of two, see eq.(8). 

The ratio Fk(Q 2 )/ F n (Q 2 ) measures the 577(3) symmetry breaking effects. 
It was obtained in [8. that the kaon wave function 4>k{x) is narrower than 
(f) n (x). This decreases the value of the integral in eq.(2) and compensates 
for fx/fn — 1-2. So, while the naive estimate looks as: F K {Q 2 ) / F^iQ 2 ) = 
(fx/U) 2 ^ 1.5, it was predicted in [12] : F K {Q 2 ) / F^Q 2 ) ~ 0.9, and this 
agrees with the recent CLEO data at s Q = 13.48 GeV 2 [T5] . see fig.2 : 

SoFk(so) = (0.84 ± 0.05 ± 0.02)GeV 2 , 

SoF^So) = (1.01 ± 0.11 ± 0.07)GeV 2 . (9) 

4. The form factor 7*771-°. 



6 The appropriate choice of a s and Jl in eq.(2) serves to diminish the role of higher 
loop corrections to the Born term. Clearly, the wider is the pion wave function ^(x,]!), 
the smaller is the mean virtuality ~p 2 ~ k 2 of the hard gluon in the lowest order Feynman 
diagram for Ft T (Q 2 ), and the larger is a s (k 2 ). 



QlFAQl) 
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Figure 2: Compilation of the existing experimental data for the pion, kaon and 
proton form factors with timelike momentum transfer. The solid points are from 
identified and K^. The open points are from unidentified h , divided into tt^ 
and according to VDM. The open triangle is from J/ip — > 7r + 7r~, supposing 
it proceeds through J/tfi — > 7* — > ir + ir~ only. The dashed curve is the fit to 
the spacelike form factor of the proton. The new CLEO points ^S] are at s = 
13.48 GeV 2 . 
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This form factor has been measured by CLEO ^B] in the range 2 GeV 2 < 
Q 2 < 8 GeV 2 . The measured value at Q 2 a = 8 GeV 2 is: 

$ 7 ,(Q 2 ) = Q 2 F^(Ql) = (16.7 ± 2.5 ± 0.4) • 10~ 2 GeV. (10) 

From the theory side, see e.g. 17 , [12J : 

§ 1 M)=h(l + 0{a s )+0{tf QCD /Q 2 )^ I =^Ml£±^{ X} ti). (11) 

For the asymptotic wave function (f)^ sy (x) = 6x(l — x) : I^ sy > = y/2 f w ~ 
0.184 GeV. This is a reference point. For Q 2 = Q 2 : a) when using (f)^ z (x, ~p, = 
0.5 GeV) in eq.(ll), it gives lf z & = (5/3) ~ 0.3 GeV; b) when using 
^ z (x,Jl ~ 2 GeV), it gives lf z & ~ (4/3)7^) ~ 0.25 GeV. It is seen that, 
unlike the charmonium decays (see eqs.(6,7)) where the decay amplitudes 
differ by a factor ~ 5 for <j)^ sy ' and <\f^ z , the difference here is only ~ 50%. 
So, to infer some conclusions about the form of 4> % {x,ii) from the CLEO 
results, it is of crucial importance to estimate reliably the loop and power 
corrections in eq.(ll). 

As for the perturbative loop corrections, the most advanced calculation 
has been performed by P. Gosdzinsky and N. Kivel JH] for (f> w (x) = (p^ asy ^(x). 
Their result (calculated in the b Q 1 - approximation) looks as: Q^ y \Q 2 ) ~ 

J^)(l-0.30) ~ 0.13 GeV. 

As for the ~ Q~ 2 power correction, the value obtained by A. Khod- 
jamirian fo r the twist-4 two- and three-particle wave functions contribu- 
tion is: A$ 17r (Q 2 ) ~ V2fn(-3 8/Q 2 ) = V2f 7T (-0.QGeV 2 /Q 2 ). It is seen 
that it has a typical value ~ ±(1 GeV 2 /Q 2 ), expected for a power correction. 
So, for 4> w (x) = ^ asy \x) one has on the whole at Q 2 = 8 GeV 2 : 

$% v) (Ql) ~ (1 -0.30 - 0.075)/^ = 0.625 -0.184 GeV = 0.115 GeV, (12) 

which looks somewhat small, see eq.(10). 

As for (f>cz(x,Ji), supposing the approach of [IE] will give the relative val- 
ues of loop corrections similar to those for (fi( asy \x), one obtains: a) when 
using <fiC z (x,~p = 0.5 GeV) for I in eq.(ll), this will give: $^ Z,1) (Q 2 ) ~ 
0.67 ■ 0.3 GeV = 0.20 GeV; b) when using ^ Z ' 2 (x,/I ~ 2 GeV), this will 
give : $^ Z ' 2) (Q 2 ) ~ 0.64 ■ 0.25 GeV = 0.16 GeV . In comparison with the ex- 
perimental result in eq.(10), all this is not worse, at least, than for (f)^ asy \x). 
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5. The new non-local axial anomaly. 



It was claimed by D. Melikhov and B. Stech in [20] that the anomaly of 
the axial current is not exhausted by the standard local terms, but contains 
a series of additional new non-local terms: 

du(rtu1*q) = 2m (qij 5 q) + Q 2 q ^(F, u F; u ) + ^.((3^G^) + 



+ eQ g k — 



1 

2 



+ 0(l/d 4 



(13) 



where the constant 

The explicit example considered in [20] from which the above eq.(13) has 
been inferred, was the form factor 



T u = (p-(q 2 )\A u (0)\ 1 (q 1 )) 



(14) 



where: A v = d^ v ^u is the axial current with the momentum p = q2 — qi, 
7(<7i) is the on-shell photon and p~{q?) is the p~ - meson. This form factor 
is very similar to the form factor 7*771° considered in the previous section. 

It has been obtained in [201 that (in the chiral limit and at large p 2 ) the 
divergence of A u is non-zero, see fig. 3: 



Vv (p (q 2 )\A v \n/( qi )) ~ K 



p2 



(15) 



and the constant k was calculated explicitly through the integral over the p 
- meson wave functions (see jT^j for the definitions and asymptotic forms of 
the p- meson wave functions, we use here the same notations and the same 
asymptotic forms of these wave functions as in [20], $(x) = 3x(l — x)(2x — 
l)/2 , V ± (x) = 3[1 + (2x - l) 2 ]/4 , V A (x) = 6x(l - x) ) : 



1 dx r $(x) - (1 - x)V ± (x) t (l + x)V A (x) 



x 



Ax 2 



(16) 



Substituting into eq.(16) the above explicit form of the wave functions, one 
obtains: k = —3/2 7^ 0. This is the main result of the paper [2*U] . 

This result looks highly surprising because the anomalous contribution 
originates here directly from the original light quark operators, see fig. 3 . 
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p(q 2 X 




p(q 2 ). 



Figure 3: The diagrams for the new non-local anomaly 

One way to make the axial anomaly "visible", is to use from the beginning 
the heavy regulator quark fields Then the equations of motion are the 
standard ones: ^^7^755+^7^75^ = 2m qi'-f 5 q + 2M ^75^/. But the heavy 
regulator fields \1/ are absent in the external low energy states. So, the reg- 
ulator fields have to be contracted into the loop, from which only the gauge 
fields can originate. This gives the standard form of the axial anomaly. So, 
the light quarks do not contribute directly to the anomaly (at vn u ^ = and 
/ 0), and the constant k in eq.(16) should be zero. On the other hand, 
the calculation of k = —3/2 from j20j l°°ks right, at the first sight. So, what 
is going wrong ? 

Our answer is that the right form of eq.(16) looks really as: 



K 



lim / dx 

s^o Jo 



$(x) - (1 - x)V±{x) (1 + x)V A {x) 



(x + 5) 



A(x + 5) 2 



(17) 



where 5 is a small power correction ~ Aq CD /p 2 , which is present in the quark 
propagator in fig. 3 (because quarks inside the p - meson are not strictly on- 
shell, but have virtualities ~ Aq CD ). Calculating k from eq.(17) one obtains: 
k = 0. 

The difference with originates clearly from the fact that in (20] the 
order of lim^o and Jq 1 dx was interchanged, but this is not allowed in this 



case. 



We conclude that there is no any new non-local axial anomaly. 



6. Large angle cross sections 77 — > n + n , K + K , K$Ks 



7 Another way, one can first integrate by 
The contribution of the total derivative is zero 
this (5 can be put zero even under the integral, 



parts the last term with Va(x) in eq.(17). 
at 5 0, and this is a crucial point. After 
and one obtains n = after integration. 
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Figure 4: Two typical Feynman diagrams for the leading term hard contributions 
to 77 — > MM , the broken line is the hard gluon exchange. 

The leading contributions to the hard kernels for these amplitudes at 
large s = W 2 = (qi + q 2 ) 2 and fixed c.m.s. angle 9 were calculated first in 
21j for symmetric meson wave functions, 4>m{ x ) — 4>m (1 — x )i an d later in 
[22] (BC in what follows) for arbitrary wave functions. Two typical Feynman 
diagrams are shown in fig.4 . The main features of these cross sections are as 
follows (below we follow mainly the definite predictions from BC in [22]) • 
a) do"(7T + 7r~) can be written as: 

s 3 da(77 - vr+TT-) _Wf\s,9)\ 2 -\ sF * ead) ^\\i- v{ e)\\ (18) 



lQna 2 d\cos9\ sin 9 sin 9 

where F^ ead ^ (s) is the leading term of the pion form factor pQ : 

\sFi lead \s)\ = \Uj o —Mx, V)\ • (19) 

While the value of F^ ead \s) is very sensitive to the form of (f) n (x), the 
factor v(9), as emphasized in [21], is nearly independent on the form of <j) n (x), 
and depends only weakly on 9. Numerically, v(9) ~ 0.12. 

The recent data from the Belle collaboration [22] for (7r + 7r~) and (K + K~) 
agree with ~ 1/ sin 4 9 dependence at W > 3 GeV, while the angular distri- 
bution is somewhat more steep at lower energies. The energy dependence at 
2.4 GeV < W < 4.1 GeV was fitted in 23, as: ^(tt+tt") = / a6 dc(da/d\c\) ~ 
W~ n , n = (7.9 ± 0.4 ± 1.5) for (tt+tt"), and n = (7.3 ± 0.3 ± 1.5) for 
(K + K~). The overall value n ~ 6 is also acceptable however, see fig. 5 . 
As for the absolute normalization, the (tt + 7t~) data are fitted in (22! to the 
eq.(18) with: \^ f ^(s,9)\ = (0.503 ± 0.007 ± 0.035) GeV 2 . 8 This value 
can be compared with : 0.88 • \sF^ cz \s)\ ~ 0.4 GeV 2 for 4>^ z \x, fi ), and 

8 Clearly, in addition to the leading terms A^ lead ' , this experimental value includes 
also all loop and power corrections 8 A to the 77 — ► 7r + 7r~ amplitudes A — A^ lead ^ + 5 A. 
These are different of course from corrections 5F W to the genuine pion form factor F w = 
p(iead) g Q ^ direct connection between the leading terms of <ier(7r + 7r~) and l-F^) 2 

in eq.(18) does not hold on account of corrections. 
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0.88 • \sF( asy \s)\ ~ 0.12 GeV 2 for ^ asy \x). It is seen that the wide pion 
wave function <p^ z \x) is preferable, while 0( asj/ ) (x) gives the cross section 
which is ~ 15 times smaller than data. It seems impossible that, at energies 
s = 10 — 15 GeV 2 , higher loop or power corrections can cure so large differ- 
ence. 9 



b) The SU(3)-symmetry breaking, do(K + K~) ^ dcr(7r + 7r~), originates 
not only from different meson couplings, fx ^ fn, but also from symmetry 
breaking effects in normalized meson wave functions, 4>k{%) 7^ 4>n{x). These 
two effects tend to cancel each other. So, instead of the naive prediction 
— Wk/ fir) A — 2.3 from [21], the prediction of BC for this ratio is close to 
unity, and this agrees with the recent data from Belle j2S| : 



<r (77 — > K + K 



Uk/U) 



2.3 



1.06 



Brodsky, Lepage 



Benayoun, Chernyak [22 



_ (0.89 ±0.04 ±0.15) Belle -2004 



c) The leading terms in cross sections for neutral particles are much 
smaller than for charged ones. For instance, it was obtained by BC that 
the ratio da ilead \7r 7c )/da {lead \n + n-) varies from ~ 0.07 at cos9 = to 
~ 0.04 at cos 9 = 0.6. Besides, it was obtained therein for the ratio: 
(K°ir o )( /ead V(7r 7r°) (W) - 1-3 • (4/25) ~ 0.21 . So, for instance, one obtains 
for cross sections <j( Zead ) integrated over < | cos#| < 0.6 for charged particles 
and over < cos# < 0.6 for neutral ones : a^ ead) (K s K s )/a^ ead) (K + K^) ~ 
0.005 . 

It is seen that the leading contribution to a Q (KsKs) is very small. This 
implies that, unlike to the case a (K + K~), it is not yet dominant at present 
energies W 2 < 16 GeV. In other words, the amplitude ^(77 — > K S K S ) = 
(a(s,9) + b(s,9)) is dominated by the non-leading term b(s,9) ~ g(9)/s 2 , 
while the formally leading term a(s, 9) ~ C Q fBc{@)/ s has so small coefficient 

9 A similar situation occurs in calculations of charmonium decays. Br(x Q — > n + ir~) 
and Br(x2 tt + tt~ -) calculated with (/> n (x) = (j) (asv) (x) are ~ 20 times smaller than the 
data, while the use of 4>tt(x) — $- cz ^ (x, /z Q ) leads to values in a reasonable agreement with 
the data, see eqs.(6, 7). 
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W[GeV] 



Figure 5: The cross sections a Q = J dcos 6 (da /dcos 9) , integrated over the cm. 
angular region |cos#| < 0.6, for a) 77 — ► 7r + 7r~ , b) 77 — ► K + K~ , together 
with a W~ 6 dependence line ; c) the cross section ratio, the solid line is the result of 
the fit for the data above 3 GeV , the errors indicated by short ticks are statistical 
only. 



15 




Figure 6: a) the overall picture of the handbag contribution, b) the lowest order 
Feynman diagram for the light cone sum rule 

C Q that \b(s,8)\ > \a(s,0)\ at, say, W 2 < 12GeV 2 . So, it has no mean- 
ing to compare the leading term prediction of BC (i.e. da(KsKs)/d cos 9 ~ 
\a(s,9)\ 2 /s ~ \ fBc(9)\ 2 /W 6 at s — > oo) for the energy and angular depen- 
dence of da(KsKs) with the recent data from Belle [21]. Really, the only 
QCD prediction for 6GeV 2 < W 2 < 12GeV 2 is the energy dependence: 
da(K s K s )/dcos9 ~ \b(s,9)\ 2 /s ~ \g(9)\ 2 /W 10 , while the angular depen- 
dence |<7(#)| 2 and the absolute normalization are unknown. This energy de- 
pendence agrees with (21j, see fig. 7. 

The hand-bag model (DKV in what follows) is a part of a gen- 
eral ideology which claims that present day energies are insufficient for the 
leading terms QCD to be the main ones. Instead, the soft nonperturbative 
contributions are supposed to dominate the amplitudes. The handbag model 
represents applications of this ideology to description of dai^yj — ► MM). It 
assumes that the above described hard contributions really dominate at very 
high energies only, while the main contributions at present energies originate 
from the fig. 6a diagram. Here, two photons interact with the same quark 
only, and these "active" gg-quarks carry nearly the whole meson momenta, 
while the additional "passive" q'q' quarks are "wee partons" which are picked 
out from the vacuum by soft non-perturbative interactions. It was obtained 
by DKV that the angular dependence of amplitudes is ~ 1/ sin 2 9 for all 
charged and neutral pions or kaons, while the energy dependence is not pre- 
dicted and is described by some soft form factors Rm(s), which are then 
fitted to the data. Because the "passive" quarks are picked out from the 
vacuum by soft non-perturbative forces, these soft form factors are power 
suppressed at sufficiently large s : Rm(s) < l/s 2 , in comparison with the 
leading meson form factors, Fm(s) ~ l/s. 

Some specific predictions of the handbag model look as : 
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n°n°)/a (n- 


-TP") 


= 1/2, and 








°o{K s K s ) 


1 


e 2 R s ^ d (s) + e 2 d R d ^ 


,00 


2 _ 2 


1 - A 


a (K+K~) 


~ 2 


e 2 u R u ^ s (s) + e 2 R s . 




25 


1-iA 

5 



> 0.08. (20) 



Here, e u = 2/3, = e s = — 1/3 are the quark charges, while the form factor 
R u ^ s (s) corresponds to the active w-quark and passive s-quark, etc. It seems 
clear that it is harder for soft interactions with the scale ~ ^qcd to pick out 
from the vacuum the heavier ss-pair, than the light uu- or dd-pairs. 10 So : 
R u ^ s (s) I R s ^ u (s) = (1 — 2A) , A > 0. (The same inequality A > follows 
from the fact that the heavier s-quark carries, on the average, the larger 
fraction, (x s ) > 0.5, of the K-meson momentum). Therefore, the handbag 
model predicts that the number 2/25 = 0.080 is the lower bound in eq.(20). 

The cross section da (KsKs)/d cos 9 has been measured recently by the 
Belle collaboration [25- The energy dependence at 2.4 GeV < W < 4.0 GeV 
was found to be : a (K s K s ) ~ W~ k , k = (9.93 ± 0.4), see fig.7 . As for the 

10 The effect due to m s ^ of the hard quark propagating between two photons in 
fig. 6 is small and can be neglected, see |2(j) . 
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angular distribution, it is consisted with da(KsKs)/dcos9 ~ l/fW 10 sin 4 9), 
|24j . The measured ratio (K $K $) / (K + K~) decreases strongly with increas- 
ing energy, and becomes smaller than the lower bound 0.080 in eq.(20) at 
W ~ 2.8 GeV. This is in contradiction with the hand-bag model predic- 
tions. 11 

Moreover, the recent explicit calculation of the hand-bag diagrams in j2H] 
using the method of the light cone sum rules [27] [2H], see fig. 6b , shows that 
for all channels, 7r + 7r~ , K + K~ and KsK$ '■ a) the energy dependence of 
the handbag form factors Rm(s) is ~ (s a /s) 2 already in the energy interval 
2.5 GeV < W = y/s < 4 GeV where the experiments have been done, and so 
this disagrees with the data on a (7r + 7r _ ) and ao(K + K~) ; b) all handbag 
amplitudes do not depend on the scattering angle 6 (in contradiction with 
the DKV results) : 

da handba9 (MM) \R M (s)\ 2 const . 

~ ~ . (21) 

dcos9 s W 10 ' K J 

and this disagrees with the data which show (~ 1/ sin 4 #) dependence; c) 
and finally, the absolute values of all three cross sections predicted by the 
hanbag model are much smaller than their experimental values. 

The conclusion for this section is that the leading term QCD predictions 
for 77 — > MM are in a reasonable agreement with the data (but only for the 
wide pion and kaon wave functions, like 4> n K (x)), while the hand-bag model 
contradicts the data in all respects : the energy dependence, the angular 
distribution and the absolute normalization. 
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